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We consider chiral fermion confinement in scalar thick branes, which are known to localize gravity, 
coupled through a Yukawa term. The conditions for the confinement and their behavior in the thin- 
wall limit are found for various different BPS branes, including double walls and branes interpolating 
between different AdSs spacetimes. We show that only one massless chiral mode is localized in all 
these walls, whenever the wall thickness is keep finite. We also show that, independently of wall's 
thickness, chiral fermionic modes cannot be localized in dSi walls embedded in a M5 spacetime. 
Finally, massive fermions in double wall spacetimes are also investigated. We find that, besides the 
massless chiral mode localization, these double walls support quasi-localized massive modes of both 
chiralities. 
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(N ' I. INTRODUCTION 

c \ 

It is a well-known fact that fermions cannot be localized in Randall-Sundrum [1] branes. The fermion chiral modes 
' ' ■ turn out to be proportional to the inverse warp factor, so that the same effect that allows to confine gravitons forces 
the fermions out of the brane. This can be avoided if the brane is in fact a domain wall generated by the vacuum 
expectation value of a scalar field, as suggested long time ago by Rubakov and Shaposhnikov [2], and fermions are 
! ' coupled to the scalar field. Bajc and Gabadadze [3] showed how this can be achieved by introducing a Yukawa coupling 
I with the scalar field, assuming an infinitely thin kink profile for the wall. 
__ ■ However, in order to have a consistent kink profile, the scalar field should be a solution to the coupled Einstein-scalar 
[ field system with a suitable symmetry breaking potential and the thin brane must be obtained as the thin wall limit 
, of this solution. Thin wall geometries have distribution-valued curvatures whose singular parts are proportional to a 
Dirac distribution supported on the surface where the wall is localized, and it is known that strong conditions must 
be imposed on a spacetime metric in order to ensure that its curvature tensors converge (in the sense of distributions) 
to the curvature tensors of the limit metric [4] . One of the relevant consequences of these restrictions is that domain 
wall solutions cannot be made infinitely thin while keeping the asymptotic values of the scalar field fixed, as in the 
step- function kink of [3] . 

Localization of gravity on thick domain walls has been considered in various works [5, 6, 7, 8, 9, 10, 11] and the thin 
wall limit (in the above mentioned distributional convergence sense) in [11, 12, 13]. Some of these thick branes reduce 
to the Randall-Sundrum (RS) brane [1] in the thin-wall limit, but not all of them, notably dynamic walls [12, 14], 
walls that interpolate between AdS§ spacetimes with different cosmological constants [11] and double- walls [10, 13]. 
In some cases, as in the asymmetric branes found in [15], localization of gravity is not possible [11], although the 
scalar field behaves as a domain wall and the stress-energy tensor is distributionally well defined even for thickness 
5.^ ■ zero. 

It is therefore of interest to investigate the possibility that fermion confinement, being directly dependent on the 
scalar field solution and not only on the spacetime metric, can be affected by the internal structure on the thick 
brane. Confinement of fermions in general spacetimes has been addressed for example in [16]. Localization of chiral 
fermions in the regularized RS scenario of [6, 7, 8] in [7]. In [17], it is found that allowing for an 5-dimcnsional 
fermion mass term can result in the confinement of both left and right modes, although a marginal one. Namely, 
fermions can escape into the bulk by tunneling, and the rate depends on the parameters of the scalar field potential. 
Confinement of massive fermions on a numerical thick brane solution, has been considered in [18], introducing an 
effective 4-dimcnsional mass term on the brane. In [19], the same results are found using the analytical solution of 
[5]. 

In this paper we address the issue of chiral fermion mode confinement in thick branes, in particular those that 
exhibit an interesting internal structure. After reviewing this confinement in two BPS regularized versions of the 
RS brane [7, 19], we consider the asymmetric BPS branes [11], where the wall separates spacetimes with different 
cosmological constants, and find the conditions for fermion mode localization. The Yukawa coupling required to 
localize chiral fermions is shown to diverge in the thin-wall limit for all these BPS branes. We then turn to spacetimes 
with a 0IS4 thick brane embedded in a M§ bulk, whose four dimensional version was found in [14]. These are known 
to localize gravity [9, 22] and could be of cosmological interest. However, we find that these walls cannot localize 
fermion modes, independently of the wall's thickness. Finally, we explore the double walls found in [13]. We show 
that massless chiral fermions can be confined in these walls, and that in the massive case they can allow for left- and 
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right-handed mode quasi-confincmcnt. 



II. LOCALIZATION OF MASSLESS FERMIONS 

Let us first consider a domain wall solution of the coupled Einstein-scalar fields equations, with an asymptotically 
AdSr, metric of the form 



ds 2 



e 2A{r) r ]iiV dx ll dx v + dr 2 



(1) 



where r\ is the 4-dimensional Minkowski metric (in the following a capital index like M runs from to 4 while 
greek indices run from to 3, we use units where 8ttGq = 1, with G5 the gravitational constant in five dimensional 
spacetime). 

The 5-dimensional spinor field satisfies the massless Dirac equation in the background (1) 



r M V M ^(x,r) = 



(2) 



With the decomposition 



*(x,r)± = ip{x) ± f{r)-_ 



(3) 



and demanding ip{x)± satisfy the massless 4-dimensional Dirac equation, the equation for the chiral modes is integrated 
to give [3] 



/(r)± oc e- 2 ^) 



(4) 



So that the very fact that the gravitational modes are confined to the wall, i.e. e 2Aij ~^ — > as |r| — > 00, implies that 
f{r)± is not normalizable, i.e. that the fcrmionic fields cannot be confined. In order to have fermions on the wall, 
one can couple them to the scalar field with a Yukawa term [3] of the form A^^, and then get 



/(r)±c<e-^» ±A ^< r '* 



(5) 



For appropriate values of A, the function /(r)_ can in principle be normalized and one chiral fermion mode confined. 
Let us consider some specific solutions for the domain wall <fi(r) 



The Kink 



The solution 



(fir \ 



(6) 



for the Einstein-scalar field system with a symmetry breaking potential V(4>) is found when [6] 



A(r) 



In ( cosh 



or \ \ I 



■ tanh" 



fir 



(7) 



with 



V(<f>) = J -a 2 



by 3 



1 - 



(8) 



and the AdS 5 cosmological constant is given by A = — 8ti 2 /3. Notice that although V((j)) as given by (8) is unbounded 
from below, a stability argument exists for this kind of BPS walls [20]. Confinement of gravitational modes for this 
model, in different parameterized forms, is discussed in [7, 8]. 

The thin-wall limit is found for 6 — > and corresponds to A(r) — > — 2a|r|/3, which behaves as the usual warp 
factor of the RS geometry, while <fi — » in this limit. Actually, following [12], it can be proved that (1,7) provides a 
sequence of metrics that satisfies the required convergence conditions of [4]. Then the distributional limit S — > of 
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all the curvature tensor fields of (1,7) exists and gives the curvatures of the limit metric. This rather technical proof 
will not be presented here. 

Now, for this domain wall, it is readily found using (5) that one chiral fermion mode is localized on the brane when 



A > 



(9) 



This result, in a slightly different parameterized form, has been found in [7]. 
there is no confinement in the infinitely- thin wall. 



It should be noted that, strictly speaking, 



B. Domain wall for a sine-Gordon potential 

A solution with 

rh = arrrtam I sinVi ( 



3i5 arctan ( sinh ( — ) ) (10) 



can be found [5] for the metric (1) with 

A(r) =-<nn(cosh(y)) (11) 



and the scalar field potential 

' 1 



V = 3a z 



2 ) cos 2 



2d ) Vv 7 ^ 



(12) 



The parameter a is given by the cosmological constant, A = —6a 2 . This smooth domain wall geometry (parameterized 
in a sligtly different form) has been shown to localize four-dimensional gravity on the wall [5] . 

The thin- wall limit for this solution is found for S — ► and corresponds to A{r) — > — a\r\, the usual warp factor 
of the RS geometry. It can be shown that in this limit the energy-momentum tensor and all the curvatures converge 
rigourously to the corresponding tensor distributions associated to the RS thin brane geometry [12, 21]), with singular 
parts which are proportional to a ^-distribution supported on the wall's plane. Now, it should be noticed that the 
scalar field (10) vanishes everywhere in the thin- wall limit S — ► 0. Since 6 is the only parameter that can be consistently 
interpreted as the wall's thickness, it is not possible to take the thin-wall limit rigorously while keeping <f> ^ 0. 

Now, equation (5) gives 

f , (ar\\** ( s-5 3 / 2 ^ E k [arctan(sinharM)] 2fc+2 1 
/W± oc(cosh(-)) cxp<j ±Av / 3 — g (2fc V 2 )(2ib)! [ 



| (13) 

where E k are the Euler's numbers (E± = 1, E<i = 5, E3 = 61, etc. and Eq = 1), which can be rewritten as 



f(r)± oc cxp < ^2 



S2 2k (2 k -l)B k ^ A^^Eh 



k(2k)\ 



± 



2ak(2k- 2)! 



[arctan(sinh ar/S)] 



2 k- 



(14) 



where B k are the Bernoulli's numbers (B\ = 1/6, B2 = 1/30, B3 = 1/42, etc.) From (14) it follows that in order to 
confine one of the chiral modes, one must require 



A > 



2a 



3 V35 



1 



2\A[ 
5 



(15) 



Clearly in this model, as in the previous one, there is no confinement in the infinitely-thin wall. 

Localization of fermions on this thick brane, with a different parametrization and employing the asymptotic values 
of 4> as r — > ±00 [3] in (5) instead of evaluating explicitly the integral as we have done, is discussed in [19]. Here it is 
instructive to compare the result obtained in this way with the exact result (15). In the large r limit, we find that 
asymptotically 



f(r)± oc exp (2a ± A V36 - I 



(16) 
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so that in order to confine one of the chiral modes, one must have 



4a 2 /2IAI , . 

Notice that solutions to the Einstein-scalar field equations with a plane-parallel symmetry with a given V(4>) and 
<j>(r) are always part of a family of "shifted" solutions with 4>(r) = <j)(r) — e and V(4>) = V(<p + e), with e a constant, 
having the same metric and cosmological constant. The shifting is irrelevant for the gravitational properties, but not 
for the fermion confinement, as is easily seen by setting e — > <p{r = oo): fermion modes can then escape to infinity, 
or in other words, the Yukawa coupling required is infinite. In the case of wall separating spacetimes with the same 
cosmological constant, the Yukawa coupling required is minimal for e = 0, the symmetric kink solution. But this is 
not the case for asymmetric walls, as we consider now. 

C. Asymmetric walls 

An example of an asymmetric domain wall was found in [11], with 

(j) = 2V35 |exp{-e- /3r/ ' 5 } - el (18) 



V{<j>) = 18 



1 - In 



125 



2V35e 



and a warp factor 



A{r) = ar - S exp(-2 e - 0r / 5 ) + <JEi (-2e~ 0r l 5 \ 



(19) 



(20) 



where the exponential integral is given by Ei(u) = -/"e^/rrfr (sec [11] for details). This is a family of domain 
wall spacetimes without reflection symmetry along the direction perpendicular to the wall. Far from the wall, the 
warp factor is 



A(r 



-00) = ar, A{r — > 00) = — ((3 — a) r 



(21) 



The wall thus interpolates between AdS$ spacetimes with cosmological constant A_ = —6a 2 for r < and A+ = 
— 6(/3 — a) 2 for r > 0. The constant e fixes the asymptotic values of the scalar field at r = ±00 and 5 is the wall's 
thickness. This asymmetric domain wall geometry (for /3 > a > 0) has been shown to localize four-dimensional gravity 
on the wall [11]. 

Now, equation (5) gives 



/(r)± oc -2 \ar - 5 exp(-2 e -' 3r/ ' 5 ) + <5Ei (-2 e - pr/5 ^ 



±2XV3S 



,-Pr/S 



(22) 



whose behavior is however rather complicated. Again, one can study the asymptotic behavior of the chiral modes. 
We have for the <fi field 

(23) 
(24) 



4>{r -> -00) = -2V3(5e <f>(r -> 00) = -2V35(1 - e) 
The Yukawa coupling required to confine fcrmions to these walls is minimized when 



e-^l + Vra/lAl 



and the condition for this confinement reads 



A > 



1 



3V2S 



(vl^M + v^M) 



(25) 



Since the scalar field depends on (3 oc ^/|A + | + ^/|A_|, fermions are confined to the wall as long as one of the 
cosmological constants is non-zero; of course gravity is not localized if A + = or A_ = 0. 
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D. A de Sitter thick brane 



A domain wall solution interpolating between spacetimes without a 5-dimcnsional cosmological constant can be 
found if one allows for a de Sitter expansion on the wall plane, i.e. a dynamic solution [14], that can be of cosmological 
interest. The Einstein-scalar field system can be integrated for a metric 



where i — 1, 2, 3. 

The scalar field is given by 



with 



dg 2 = e 2A(r) (_ dt 2 + e 2fH fai fax + df 2 



(r) = y/36(l — S) arctan f sinh 



A(f) 



S In cosh 



Qf 



(26) 



(27) 



(28) 



and a sine-Gordon scalar field potential. These walls have a well-defined distributional thin limit S — > 0, a thin dS± 
brane embedded in a M5 bulk [12, 21]). Although they do not become an RS brane in the thin wall limit, they can 
localize gravity with a mass gap between the zero modes and the massive ones [9, 11, 22]. However, they cannot 
confine fcrmions as we now show. 

In these coordinates, the Dirac equation results in a different expression for f± 



f± ocexp j-2A(f)±A J , 



<j)(r)dr 



(29) 



which can be explicitly integrated and yields a power series of arctan(sinhar/(5) as in (14), but quite more involved 
with coefficients which are polynomials in 5. Each term of the series, for increasing power, pushes up the lower bound 
on A necessary to obtain the desired normalizability of f±. Hence, there is no value of A that renders f± normalizable. 
This is not at all unexpected, in this model A(f) is the same function A(r) of the domain wall with a sine-Gordon 
potential and f± can be rewritten as 



f± oc exp J drT(r)± 



where 



Y( r )± = 2 a tanh — ± A y/3S(l — S) exp(— <51n(cosh ^-)) arctan(sinh — ) 
Id 



(30) 



(31) 



It follows that the Yukawa coupling contributes to T(f)± with a term that becomes exponentially suppressed with 
respect to the one of the domain wall of the sine-Gordon potential, becoming unable to cancel the coming one from 
cxp{— 2A(f)} for \r\ — > 00. Asymptotically, 



T(f)_ - ± (2a - Xy/3S(l-S)e-^\ | 



(32) 



as r — > ±00 and the term coming from the Yukawa coupling goes to zero. This behavior has also been checked by 
numerical integration of (29). 



E. Double walls 



In [13], a two-parameter family of walls are found with a metric 

ds 2 = e 2A(f) (rj^dx^dx" + dr 2 ) 

The scalar field solution is 

</>(f) = 4>o arctan(af) s ; 0o = ~ 



(33) 



(34) 
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with 



and 



A(f) = --ln(l + (af) 2 *) 



V = 3a 2 sin(0/0 o ) 2 - 2 / s 



2s + 3 



cos 2 (</>/</> )-2 



(35) 



(36) 



If s is an odd integer, the field (34) interpolates between two global minima of the potential, but for s > 1 there is a 
local minima between any two global ones. The configuration then is that of a double wall, with the energy density 
peaking at two points, interpolating between AdSs asymptotic vacua with A = —6a 2 . These walls have been shown 
to confine gravity in [11]; similar double walls have been found in [10]. For s = 1, this is just the RS regularized 
spacetime considered in II B, for a fixed thickness 6=1 and written in conformal coordinates. 

The chiral modes take again the form (29), but in contrast with the dynamic walls, confinement of chiral fermions is 
possible in this case, for any odd value of s. The integration of (29) in this case has again been performed numerically, 
giving a rcnormalizable mode. But it is also possible in this case to examine the asymptotic behavior of the integrand, 
writing the modes as in (30) we find 



T(f) : 



2 7T s /6s~- 
— zb A— 

f 2 s 



3 1 



ar 



(37) 



as Irl — > oo. So that the Yukawa term dominates when 



A > 



4a s 
7T V3 V2s - 1 



(38) 



III. LOCALIZATION OF MASSIVE FERMIONS 



Up to now, we have supposed that the 4-dimensional spinor fields satisfy the massless Dirac equation on the wall. 
The inclusion of a mass term that connects left and right-handed modes could modify the well-known result of having 
only one chiral mode localized on the branc. This possibility has been explored in Rcf. [18], using a numerical solution 
for the domain wall. If instead of this one uses the exact solutions (10,11) as in Ref. [19], or (6,7), very similar results 
are found. Namely, the chiral modes can be shown to satisfy a Schrodingcr-like equation, with different potentials for 
left and right-modes. These potentials are unbounded from below, but have minima where the modes can be confined, 
and the tunneling can in principle be sufficiently suppressed by adjusting the parameters. 

Here we study the case of the double wall spacetime of Ref. [13], which includes the regularized version of the RS 
geometry considered in [5] as a particular case. These branes have the advantage of being written in a coordinate 
system that leads to a Schrodinger equation for the chiral modes with the four-dimensional fcrmion masses as its 
eigenvalues, which enable us to analyze its solutions in a simpler way. 

Consider again a spinor "J in 5 dimensions, 



ty(x,r) = ip+(x)u+(r) + ip—(x)u—(r) 
where now the 4-dimensional spinors satisfy 

Then the chiral modes follow the equations (prime denotes derivative respect to f) 



d r + 2A'(f) ± A0(f)e A(F) 



± m u± 



Setting 

one can write the Schrodinger equations 



u(f) = u(f)e- 2A(F) 



-9 2 r+V± M 



u± = m 2 u± 



(39) 
(40) 

(41) 

(42) 

(43) 



7 



r\i- 





FIG. 1: Vq M (continuous line) and Vq M (dashed line) for the double walls with s = 1 , s = 3 and s = 5 



where 

V£ M = {H(r)e A{f) f ± W(f)e^) (44) 

In figure 1 we plot Vq M for different values of s. In all cases a = 1, A = 1. 
Now, notice that (43) can be written as 

Q + Qii-=m 2 u-, QQ + u+ = Tn 2 u+ (45) 

where Q = df + A</>(f) expA(f) and Q + = —df + A 0(f) expA(f), and (43) can be recast as a SUSY quantum 
mechanics problem as in [23]. It follows that the eigenvalues of the u- and m+ modes always come in pairs, except 
possibly for the massless modes. Indeed, this pairing of the mass eigenmodes is required for the existence of massive 
fermions satisfying (40). 

As a simple check of consistency, let us consider the massless modes again. For m = 0, from (45) we find 

uq± oc exp{±A J df(f>e A } (46) 
and from (42) the result (29) is recovered. Since for 

A > 4a S (47) 

> W3 % /2T — T 

one and only one chiral mode is normalizable, this chiral normalizable mode is just the one that was considered in 
the previous section and that gives a chiral four dimensional fermion localized on the brane. 

Now let us consider the massive modes. We recall that for s = 1, this is just the second RS regularized spacetime 
considered in the previous section, for a fixed thickness 5 — 1 and written in conformal coordinates. For s = 1, we 
see from Fig. 1 that right- and left-handed modes are trapped in a metastable state with a non-zero probability of 
tunneling into the bulk. This tunneling probability depends on the value of the Yukawa coupling and could be made 
small by increasing A. For the range of parameters shown, this tunneling probability is very high, to the point that 
one can exclude localization. However, for the same set of parameters but greater values of s, where double walls 
appear, we see from Fig. 1 that Vq M develops a deep minimum in the region between the walls. Now, right-handed 
modes of small mass can be trapped here in a metastable state like the paired massive left-handed ones. As before, 
the tunneling probability can be made small by increasing A. On the other hand, increasing s has the effect of making 
each of the walls thinner, accordingly, Vq M has also thinner "walls" , which would increase the probability of fermions 
escaping. A precise calculation of the tunneling probability for these massive left and right-handed modes is in order 
here, but will only make sense in a definite model in which the parameters of the potential can be related with the 
age of the Universe. This is beyond the scope of this paper and will be dealt with elsewhere [24]. It is however clear 
that the internal structure of the wall, i.e. the appeareance of double walls, can enhance the probability of localizing 
massive fermions. 
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IV. SUMMARY 



We have studied the localization of fcrmions on thick branes with a non-trivial internal structure. We find that 
double walls can confine fcrmions of a given chirality, as is also the case with the known regularized RS branes, provided 
they are coupled to the scalar field through a Yukawa coupling that is inversely proportional to the wall's thickness. 
Asymmetric BPS walls interpolating between spacetimes with different cosmological constants can confine fermions 
in similar way, with a Yukawa coupling proportional to the sum of the square roots of the cosmological constants. In 
all these cases, the Yukawa coupling required diverges in the thin-wall limit. On the other hand, dynamic walls with 
a cosmological dc Sitter expansion are unable to localize even one chiral mode. We have considered here only the five 
dimensional version of the well-known solution of Goetz [14] , but these results are expected to hold for other dynamic 
walls such as the ones found recently in [25] ; work in this direction is now in progress [24] . 

Finally, we have also addressed the issue of confinement of massive fermions in double walls. We find that, besides 
the massless chiral mode localization, massive modes of both chiralitics can be quasi-localized on double walls, being 
"trapped" in between the walls. The tunneling probability can in principle be made as small as desired by increasing 
the Yukawa coupling and fermion masses, but for a given set of parameters it is smaller in the double walls than in the 
single ones. Thus, fermion localization is shown to be affected by the internal structure of the domain wall solutions. 
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